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ABSTRACT 



Context. The fractal structure of the interstellar medium suggests that the interaction of UV radiation with the ISM as 
described in the context of photon-dominated regions (PDR) dominates most of the physical and chemical conditions, 
and hence the far-infrared and submm emission from the ISM in the Milky Way. 

Aims. We investigate to what extent the Galactic FIR line emission of the important species CO, C, C^, and O, 
as observed by the Cosmic Background Explorer (COBE) satellite can be modeled in the framework of a clumpy, 
UV-penetrated cloud scenario. 

Methods. The far-infrared line emission of the Milky Way is modeled as the emission from an ensemble of clumps with 
a power law clump mass spectrum and mass-size relation with power-law indices consistent with the observed ISM 
structure. The individual clump line intensities are calculated using the KOSMA-r PDR-model for spherical clumps. 
The model parameters for the cylindrically symmetric Galactic distribution of the mass density and volume filling factor 
are determined by the observed radial distributions. A constant FUV intensity, in which the clumps are embedded, is 
assumed. 

Results. We show that this scenario can explain, without any further assumptions and within a factor of about 2, the 
absolute FIR-line intensities and their distribution with Galactic longitude as observed by COBE. 

Key words. ISM: clouds - ISM: structure - Galaxy: disk - Infrared: galaxies - Infrared: ISM - Submillimeter 



1. Introduction 

The importance of the interaction between the interstellar 
UV-radiation and dense clouds in the ISM, determining the 
physical and chemical conditions in the surface regions of 
molecular clouds, which are modelled as so-called photon- 
dominated regions (PDRs), has been recognized since the 
first observati ons of the dominant cooling lines in the FIR , 
[Cii] 158 urn. (iRussell at al.lll980lll9 8ll:IStutzki et al 
Mizutani et all I1994D . [Oi] 63 ^m (jMelnick et all 



in the Milky W ay, but also in star f orming regions in ex- 



ternal galaxies (jKramer et al.l 1200 5: Cont ursi et al 



Malhotra et"aIll200lHUnger et al.ll2000.: .Nikola et al.l 



2002 



1998 



Madden et al.lll997HSchilke et al.Hl993f ) . An overview of the 



1988; 



1978; 



Stacev et al.l 19831) and the mid- and high- J CO lines 
(IStorev et all Il981t IWatson et all [19851: iJaffe et all [l987t 



iBoreiko k Betzlll99l[ ). Starting from the first PDR-models, 
tailored to explain the FIR lin e emission from massive star 
forming regions such as Orio n ([Tielens fc Hollenbach|[l985t 
ISternberg fc Dalgarnall989f ) the inclusion of more details of 
the heating and cooling mechanisms as well as of the chem- 
ical network, nowadays allows the modeling over a wide 
range of physical parameters and has resulted in success- 
ful modeling of many detailed aspects of observed photon - 
dominated regions (review by ,Hollen bach fc Tielenill999f) . 
PDRs thus have proven to be a very useful concept in un- 
derstanding the mutual interaction between star formation 
and the structure of the ISM, both in individual regions 
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different models is given by the recent comparison study of 
PDR codes (R61hg et al. 2007,). 

The importance of the PDR scenario is related to the 
fact that the ISM is fractal and thus most of the mate- 
rial is close to surfaces and hence affected by UV radia- 
tion, i.e. is located in PDRs. This was realized early on by 
comparing the spatial distribution of o bserved PDR tracers 
(jStutzki et al.lll988HHo"we et al.lll99lf ) with simple models 
of homogeneous or clumpy cloud structures. The observed 
relatively uniform line ratios of low-J ^^CO and ^^CO in 
molecular clouds, inconsistent with simple, uniform cloud 
models, are shown to be naturally explained if the emission 
i s assumed to origi nate in many, relatively small clumps 
(IStorzer et al.|[l996l ). The necessity of high densities in or- 
der to explain the observed large ^^CO brightnesses in the 
mid- J line s in the submm independently indica te a clumpy 
structure (jGraf et al.lli990HWolfire et al.lll989[ ). 

The fractal structure of the ISM is well represented 
by an ensemble of clumps with a power law clump 
mass d istribution, and a power law clump mass-size re- 
lation (jStutzki et al.|[l998[ ). Power law mass spectra have 
been derived by decomposing the observed 3D-datacubes 
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of emission of CO isotopolosues into clumps by vari 
ous methods fStutzki fc Gust en|[T990l: iKramer et ai.lll998l: 



[Williams k Blitz 1993; Will iams et al 



19941 ). The identi- 



fied clumps typically show also a power law mass-size re- 
lation, although obtaining significant coverage over more 
then 1.5 to 2 orders of magnitude in length scale is difficult; 
combining low and high angu lar resolutio n observations of 
the Polaris Flare, [Heithauscn et al.l (jl998t ) demonstrated a 
power law mass-size relation over 3 orders of magnitude 
in length scale. Calculating the PDR emission from an en- 
semble of clumps with a given mass spectrum and mass-size 
relation thus offers a convenient way to model the submil- 
limeter and FIR line emission of the ISM as an UV pene- 
trated clump ensemble with the given fractal characterstics. 

The far-infrared absolute spectrophotometer (FIRAS) 
on the COBE satellite conducted a spectral line survey of 
the Milky Way in the far-infrared region a,t wavelengths 
longw ard of 100 urn. (| Wright et all 1 19911 : iBennett et all 
Il994[ ). The spatial and spectral resolutions of respectively 
7° and 0.45 cm^^ allowed detection of the CO rotational 
transitions from J = 1-0 to 8-7, and the fine structure tran- 
sitions of [Cli] 158 Atm, [Nil] 122 and 205 fim, [Oi] 146 fim, 
and [C i] 370 and 609 fim, in addition to the dust FIR con- 
tinuum. The wavelength coverage of COBE FIRAS did not 
allow to detect a few other important cooling lines of the 
ISM like the [ O i] li ne at 63 /xm or the [O ill] line at 88 (im. 
iFixsen et al.l (jl999f ) report the distribution of spectral line 
emission along the Galactic plane. They re-binned the data 
to a 5°-grid in Galactic longitude and assumed an latitudi- 
nal extension of the FIR line emission of 1° according to the 
results of the DIRBE FIR continuum observations. In a first 
attempt to interpret the observed emission with PDR mod- 
els, they fit the plane-parallel PDR model of HoUcnbach 
(1991) to the observed line ratios in the Galactic center re- 
gion to derive a densi t y of 3 cm^"^ and FUV field of 10 
Habing unitfl iHeilej (|l994[) used the COBE Milky Way 
data to analyze the prop erties of the extended ionized gas 
traced by the [N ii] lines. iMisiriotis et al.l (j2006l ) combined 
COBE DIRBE and FIRAS continuum data at wavelengths 
between 1 /xm and 1 mm wavelength to constrain exponen- 
tial axisymmetric models for the spatial distribution of the 
dust, the stars, and the gas in the Milky Way. 

In this paper we investigate to what degree the global 
FIR-line emission distribution of the Milky Way as observed 
by COBE can be explained in the framework of an UV- 
penetrated, clumpy cloud scenario of the ISM. We mention 
that the [Ci] 1-0 and 2-1 and CO 4-3 and 7-6 emission 
from a massive star forming region in the Carina nebula can 
be explaii ied by the emission of clumpy PDRs as recently 
shown by iKramer et al.l (j2008( ). The paper is organized as 
follows: In Section 2 we shortly summarize the concepts to 
quantify the clumpy, fractal cloud structure and introduce 
the clumpy cloud PDR model. Section 3 discusses the line 
emissivity of the clump ensemble, in particular the depen- 
dence on varying mass limits of the clump distribution. In 
Section 4 the model results are compared with the obser- 
vational data from COBE. A short conclusion with outlook 
is given in Section 5. 




1 Habing unit corresponds to the integrated flux in the wave- 
length range f rom 91.2 to 111.0 nm of 1.6 x 10^'* erg s^^ cm^^ 
l|HabinElll96l ). 



-4 -3 -2 

og( Impact parameter p [pc] ) 

Fig. 1. Brightness temperature vs. impact parameter p of 
the CO 7-6 rotational transition for selected clumps with 
mass and density values following the clump-mass spectrum 
and mass-size relation, is the gas density at the surface 
of the clumps in [cm^'^] (see Eq. [T]), Af^j the clump mass in 
[Mq], and the incident FUV-flux is x = 10^ Xd- 



2. The galaxy model: an UV-illuminated clump 
ensemble 

Modeling the large-scale Galactic FIR line emission as ob- 
served by COBE is done in the following steps: 

— generate a large parameter grid for the emission of in- 
dividual spherical clumps with the KOSMA-t PDR- 
model, where the individual clumps are characterized 
by their mass, density (resp. size), and the UV- intensity 
that they are embedded in. 

— generate a clump mass and size distribution at each 
Galactocentric radius from the observed mass density 
and volume filling factor distribution in the Milky Way, 
assuming an universal power law index for the clump 
mass and size spectra and applying adequate high and 
low mass cutoffs, 

— estimate the Galactic radial distribution of the effective 
mean UV-field in which the dense ISM clumps are em- 
bedded and the average clump density in the ensemble 
respectively the volume filling factor of the clump en- 
semble, 

— by combining the above, calculate the volume emissivity 
of the ISM originating from the UV illuminated clump 
ensemble as a function of Galactocentric radius in the 
Milky Way 

— calculate the line emission versus Galactic longitude by 
line-of-sight integration over this emissivity in the di- 
rection towards the Earth. 

We describe each of these steps in detail in the following: 

2.1. The KOSMA-T PDR model: line intensities of an 
individual clump 

The KOSMA-r PDR-model describes spherical clumps 
characterized by their mass, density (or size, respectively), 
and the incident UV field. 

We treat the mass of the clumps Mc\ as the total hydro- 
gen gas mass. To obtain the total gas mass of a clump the 
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Table 1. Properties of the KOSMA-r PDR- model elumps with given masses and densities at an incident FUV-flux of 
10^ Xd- Orders of magnitude are given in parentheses. 





Densities 


Mass 


Column densities 


Temperatures 


Surface 


Center 






Clump average 


log(ns/cm 


log(nc/cm"^) 


log(Mel/M0) 




(A'co) 


(Th,) (Tco) 


3.5 


4.55 


3.0 


1.41(22) 


1.71(18) 


21.7 12.3 


4.0 


5.05 


1.5 


9.74(21) 


1.50(18) 


23.8 11.0 


4.5 


5.55 


0.0 


6.69(21) 


1.07(18) 


26.5 12.5 


5.0 


6.05 


-1.5 


4.57(21) 


7.21(17) 


29.6 15.5 


5.5 


6.55 


-3.0 


3.12(21) 


4.94(17) 


32.8 19.9 



mass of helium and heavier ele ments has to be added. Thi s 
is a factor of about 1.4 higher () Anders fc Grevessill989f ). 
Due to the fact, that the massive clumps dominate the en- 
semble mass (see Section 12. 2[) and contain only a few per- 
cent of their gas mass in the atomic phase, we can neglect 
the contribution of the atomic gas so that the ensemble 
mass is basically given by the mass of molecular hydrogen. 

The radial density structure n{r) of the model clumps 
with radius Rd is given by an inner core region with con- 
stant density for r < 0.2 Rc\ and a power law with index 
-1.5 for 0.2i?ci <r < Rei- 



ner) = 




-1.5 



, for 0.2 Rci<r < R^i 
, for r < 0.2 R^ 



(1) 



Note that rig = nn.s + 2nH2.s is the total hydrogen 
nucleus number density at the clump surface. By neglecting 
the presence of helium and metals in the number density, we 
underestimate the number of collision partners by a factor 
of 1.2. The effect on the collision rates is, however, smaller 
due to the larger mass and thus lower velocity of helium. 
This does not change the results significantly. With the 
clump mass Md, the clump volume Vc\ — ^ R^^, and the 
mass of the hydrogen atom mn , the average clump density 
Uci = Mci/(mHT4i), used in the following to characterize 
each clump, is given by rid = (2 — 5"^''^) Ug ~ 1.91 Ug. 

The FUV flux illuminating each clump is assumed to 
be isotropic. It is given in units of the Draine field xd = 
2.7 X 10"'^ erg cm~^. This is the flux of the local inter- 
stellar radiation field (ISRF) int egrated over t he wavelength 
interval from 91.2 to 111.0 nm ()Drainelll978[) . 

In the KOSMA-T PDR-model a constant flux of cos- 
mic rays and the incident FUV radiation are the only heat- 
ing sources of the clumps. The temperature and chemical 
structure of the clumps are derived from a radiative transfer 
computation including line-shielding for the incident FUV 
radiation and an escape probability approximation for the 
emitted FIR line radiation. The efficiency of the line cool- 
ing depends on the average line width Sv, which is assumed 
to be 1.2 km s~^ (see Table El). More detail s are described 
bvlGierens et all (119921) iKdster etall (|l994[ ). IStorzer et all 
(|1996D . and lRolhg et all (|2006l) ! 



In order to calculate the average line intensity Id of 
each clump, first a line-of-sight integration of the source 
function and the attenuation along parallel paths through 
the clump is performed to calculate the resulting intensities 
for different impact parameters p. As an example the 
line integrated CO 7-6 intensities are shown as a function of 
impact parameter in Fig.[T] Note the strong emission of the 
small, low mass, high density clumps. The physical reason 



Table 2. The parameter values of the precomputed 
KOSMA-r PDR-model clumps are equidistant on logarith- 
mic scale in steps of half an order of magnitude. 



Quantity 


Parameter 


range 


Density ns/cm^"* 


lo^ 


10" 


Mass Mci/M0 


lO"^ ... 


10^ 


FUV flux x/XD 


10",..., 


10" 



for this is the enhanced population of the J = 7 rotational 
excited state (i.e. a higher excitation temperature) of the 
CO molecule in the smaller and denser clumps, as only for 
those the density in the clump center falls near or above the 
critical den sity of the CO 7-6 rotational transition of about 
10^ cm-3 ([Kaufman et al.lll999t [Kramer et al.ll2004[ ). The 
less dense but massive clumps have much larger column 
densities, partially compensating for the lower excitation 
temperature and increasing the optical depth. This leads 
to the apparently constant central brightness temperatures 
of the bigger clumps in Fig. [TJ Detailed information about 
these clumps is compiled in Table [1] 

Second, Id is obtained by averaging over the projected 
clump surface via 



2 
R 



■cl Jo 



(2) 



The results of these computations are used for the clump 
ensemble intensities. 

We generate a large grid of intensities for the three inde- 
pendent clump parameters mass, density and FUV inten- 
sity, with grid points equidistant on a logarithmic scale in 
steps of half an order of magnitude. The range covered by 
the parameter cube of the KOSMA-t PDR-model is shown 
in Table [5| For clumps with intermediate parameter values, 
the KOSMA-T PDR-model intensities are interpolated be- 
tween the grid points using logarithmic differences. Further 
parameters for the PDR model, held fixed for the present 
work, are listed in Table [3] 

2.2. Clump mass distribution and geometry 

Observations of the ISM show structure at a U scales down 
to the resolution limit, a s stu d ied by e.g. Stutzki et al.[ 
(|198^,IStutzki fc Gustel] (ll990D.[Heithausen et all ([1998f ). 
[Kramer et all l[1998l) . ISimon et all ([200lh . Quantitative 
analysis of the observed intensity distribution shows that 
it can be described as a self-similar, i.e. fractal structure. 
Alternatively, the observed 3d spectral line data cubes can 
be decomposed into clump ensembles. The "clump-mass 
spectrum" , i.e. the number of clumps dA^ci per mass bin 
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Table 3. The fixed parameter values of the prccomputed 
KOSMA-r PDR-model clumps. Note that powers of ten 
are given in parentheses and the elemental abundances are 
given in units of the hydrogen abundance, i.e. Xi = rii/n. 



Quantity 



Value 



He abundance Xho 




0.1 


abundance Xo 


3.0 ( 


-4) 


C abundance Xc 


1.4(- 


-4) 


S abundance Xs 


2.8 (■ 


-5) 


^^C abundance Xisq 


2.1 ( 


-6) 


Cosmic ray ionization rate XCR [s^^] 


5(- 


17) 


Velocity width (FWHM) Sv [km s"^] 




1.2 


Dust UV cross section adust [cm^] 


1.9 (- 


21) 



dMci, is then given by 



dTVd 
dMr. 



= A An 



(3) 



The size of the clumps is correlated with their mass via the 
"mass-size relation" : 



Mr: 



CR2 



(4) 



The indices a and 7 of the clump ensemble are related 
to the fractional Brownian-motion power la w index char- 
acter izing the fractal structure of the ISM ([Stutzki et al.l 
[TUl). We choose v alues of 1.8 and 2.3, resp ectively, follow- 
ing the results from lHeithausen et al.l (|1998[ ). The ensemble 
mass Mens is obtained by integrating the clump-mass spec- 
trum over a mass range from a low mass cutoff mi to a 
high mass cutoff rriu- The total mass of the clump ensem- 
ble Mens = Mc\ jIj^ dMci determines the constant of 
proportionality A in Eq.jsjvia 



Mr- 



A 



(5) 



Note that for 1 < a < 2 the bulk of the ensemble 
mass is contained in the high mass clumps. In contrast, 
the total number of clumps A^ens = 3^ dA/ci — 

{m\^°' — mj^") is dominated by the low mass clumps. 

Using the volume Vc\ = ^R^\ of a clump of mass M^x and 
the mass-size relation, the volume of the clump ensemble 
^ens = C: 3^7 ^ci dMei is givcu by 



14, 



A 



A 1+3/7- 

47r mu 



1+3/7- 
'1 



C3/7 3 



1 + 3/7- a 



(6) 



For the values of a and 7 given above, the exponent is 1 + 
3/7 — a ~ 0.5, and hence the ensemble volume is dominated 
by the large and most massive clumps. 

The above applies as long as we are concerned with the 
integral properties of the clump distribution. In order to 
take the spatial distribution of the clumpy medium into ac- 
count, we now switch to a volume specific notation where 
the medium is characterized by the (volume) number den- 
sity, rather than the number, of the clumps, and by the 
mass density rather than the total mass. We obtain for 
the volume- and mass-differential clump distribution, i.e. 
for the (volume) number density of clumps dA'ci per clump 
mass interval dMd, the "clump number density spectrum" 
as a function of Galactocentric radius r: 



d^ci(r-) 
dMn 



i(r)A/ 



— a 
:1 • 



(7) 



The "clump mass density spectrum" , i.e. the mass per vol- 
ume element, per mass interval of the clump ensemble is 
correspondingly given by 



diVci(r) 
dA/.i 



i(r)Af, 



l-a 
cl ' 



(8) 



and the mass differential clump volume per volume element, 
i.e. the "clump filling factor spectrum" , is 



d^ci(r-) 



Ki = i(r-)C-3/T(r)A/; 



(9) 



Integration over the clump mass leads to the "volume mass 
density" p{r) and the "volume fiUing factor" /v(r), respec- 
tively, given by 



p{r) 
Mr) 



A{r) 

2- a 



47r A{r) mu 



1+3/7-0 



and 

1+3/7-a 
mj ' 



3 C^h(r) 



1 + 3/7 



(10) 
(11) 



The "ensemble averaged clump mass density" pcns{r) 
Pir)/fv{r) is given by 



sir) 



£^3/7^^Ni±^/l^ 



m 



(mu 



1+3/7- 



m 



1+3/7-Q 
1 



(12) 



As discussed in detail in subsection 2.4 below, the distri- 
bution of the dense ISM in the Milky Way in our model is 
specified by the volume mass density p{r) and the (higher) 
ensemble averaged clump mass density Pons('") (respectively 
the volume filling factor of the molecular material). For 
given values of the low and high mass cutoffs of the clump 
distribution, equations ^TU\i and p2p then determine the 
pre-factor of the clump number density spectrum and its 
spatial variation. 



A(r) = 



2- a 



rriu " — m? 



-Pir), 



(13) 



through the volume mass density profile of the Milky Way, 
and the pre-factor of the mass-size relation and its spatial 
variation. 



Cir) 



An 



2- a 



1+3/7-Q 



1+3/7-Q 
*1 



3 1 + 3/7 -a 



\ 7/3 
■Pcns(r) (14) 



through the ensemble averaged clump mass density. 

A further quantity of interest is the beam filling factor 
/b, as the ratio of the solid angle filled by the clumps of 
the ensemble Jlens and the solid angle of the beam ■ The 
solid angle of a specific clump with mass Md and Radius 
Rci at a distance s = |s| (considering s as given in Galactic 
coordinates centered at the sun) is given by 



n^iis) = 2^(1 - Vi - (Rci/s)^) 



R. 



(15) 



M. Cubick et al.: A clumpy-cloud PDR model of the global far-infrared line emission of the Milky Way 



5 



for i?ci ^ s. Then, the "solid angle density spectrum" , i.e. 
the solid angle per volume and clump mass clement of the 
ensemble, is 



dAL 



- 2/7- 



(16) 



For the adopted values of a and 7 the exponent 2/7 — 
a has a value of —0.93, implying that on a logarithmic 
mass scale, [Jci(s) dMci/d(ln Mci) = MdUJciis) oc 
each logarithmic mass interval contributes about equal to 
the solid angle of the ensemble. Integration over the given 
clump mass range gives the "solid angle density" of the 
ensemble 



,{s) = / Wci(s) dAfci 

J nil 

7/ N 2/7-Q + l 2/7-Q+l 

A{s) TTiu —my 



C2/7(s)s2 



rl/3 



2/3 



2/7 - a +1 

p{s) (2-a)i/3(i + 3/^_a)2/3 



3)2/' 



3 „2 



2/7 - a + 1 



2/7-a+i 



2/7-Q+i 
^1 



{ml 



')l/3(TOi+3/7-" _ ,„J+3/-r-«)2/3- 



(17) 



For the values of a and 7 adopted above, the nomi- 
nator shows again a very weak mass dependence of oc 
j^o.o7^ In the denominator in both factors the upper- 
mass cut-off dominates, so that the total dependence of 
the solid angle volume density of the clump ensemble is 



oc TOu 



(4/3-Q-1-2/7) 



-0.4 



i.e. weakly decreasing with in- 
creasing upper mass cut-off and basically independent on 
the lower mass cut-off. Integration over the beam volume 
Vb gives the total solid angle filled by the clump ensemble 



2/3 



(2-a)i/3(i + 3/^_a)2/3 

2/7-0:4-1 

2/7-Q+i 2/7-Q+i 

mu — TTij 

' {ml-" - m^")i/3(„ji+3/7-" - ^i+3/T-")2/3 



(18) 



with the same dependence on upper and lower mass cut-off 
as for the volume specific solid angle. 

2.3. Volume emissivity and intensity of ttie dumpy medium 

In order to calculate the intensity of the clump ensemble, 
we neglect any absorption (in particular self-absorption in 
the spectral lines) among different clumps along the line- 
of-sight. This assumption may be violated locally in the 
high densities of individual star forming regions, although 
even there, the higher velocity dispersion of the virialized 
clump ensemble helps to avoid line-of-sight crowding in each 
velocity interval. It is reasonable for the large scale FIR 
line emission of the Galaxy, taking the additional velocity 
spread between individual regions due to the differential 
Galactic rotation into account. The intensity then simply 



is given by the linc-of-sight integration over the volume 
emissivity of the clump ensemble 



/ = y 77(5) ds 



The beam average intensity then is 
^B==7r/ fri{s)dsdn = l- 



Vb 



dV 



(19) 



(20) 



The KOSMA-T PDR-model gives the clump average spe- 
cific intensity I^i for each clump of a given mass. The "vol- 
ume emissivity spectrum" then is given by multiplication 
with the clump projected area Ac\ ~ Q,ci and the clump 
number density spectrum, namely 



d77(5) _ diVel 2 



dA€ 



S Id ■ 



(21) 



With Lo, 



tJA/ci 



r^ci. we obtain for the contribution to the 



ensemble intensity in the mass bin dA/ci 
d/B 1 



dA^c 



B JVb 



/clt^cl(s) dV. 



(22) 



Integration over the mass ensemble results in the beam av- 
eraged intensity of the ensemble 



B J Vb J inj 



/cl CJcl(s)dA/cl dV. 



(23) 



Particularly, we obtain for a constant clump intensity the 
resulting beam average intensity /b with the beam filling 
factor /b = r^ens/f^B as Ib = /b ^1- 

2.4. Galactic parameter distributions 

The Galactic parameter distributions of mass, density, and 
UV-intensity are assumed to vary only in Galactocentric 
radial direction within a cylindrically symmetric Galactic 
disk of constant half-height h. The values of the different 
parameters (at the solar circle at R = R0 = 8.5 kpc) are 
compiled in Table S) We adopt a half-height of the molec- 
ular disk of 59 pc a nd for the atomic gas layer of 115 pc 
(jWolfire et al.1120031 ). The fractal characteristics (a, 7, wi, 
mu) of the clump distribution are assumed to be constant 
over the whole Galaxy. Due to the simplified structure of 
the Galactic parameter distributions in our model disk we 
do not expect to reproduce individual emission features 
along the Galactic plane. The Galactic mass distribution of 
the neutral ISM is a key parameter of the model, because 
the modeled line intensities are basically proportional to the 
mass of the molecular gas assumed to be completely repre- 
sented by PDRs (see Section [ij. The Galactic m ass distri- 
bution has been discussed in the litera t ure e.g. by|Clemensl 
Tgs l . lRohlfs k K reitschmannI p987[) . I Wilhams fc McKed 
19971 ). [Bronfman et al. (2000). We use an approximation 



of the radial Galac tic II2 mass distribution presented by 
Wolfire et all (I2003D for 3 < i?G/kpc < 18. According to 
WiUiams fc McKed (Il997l Fig. 3) the H2 mass surface den- 
sity decreases inwards from 3 to 1.7 kpc Galactocentric 
radius exponentially by a factor of about 3/8. For 0.6 
< i?G/kpc < 1.7 the molecular hydrogen gas mass distribu- 
tion is set constant. In the Galactic center region {Rq < 0.6 
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Table 4. Parameters of the Galaxy model (values at the solar circle). - * For distribution over Galactocentric r adius 
see Fig. [2] or text. - * * Assumed to be consta nt over the whole Galactic disk (see text). - [1] IWolfire et al.l (|2003f ) - [2] 
iRevnolds. (,1991. ) - [3] iHeithausen et all (|l998t ) 



Name 


Symbol 


Unit 




Value 


Reference 


Radial extension of Galactic disk 


Rg 


kpc 




18 


[IJ 


Solar distance to Galactic center 


R0 


kpc 




8.5 


[1] 


Half height H2 ** 




pc 




59 


[1] 


Half height CNM/WNM ** 


hni 


pc 




115 


[1] 


Half height WIM ** 


hwiM 


pc 




10^ 


[2] 


H2 mass surface density * 




Mq pc- 


- 2 


1.4 


[1] 


CNM mass surface density * 


ScNM 


Mq pc- 


- 2 


2.25 


[1] 


WNM mass surface density * 


SwNM 


Mq pc- 


-2 


2.75 


[1] 


WIM mass surface density * 


Swim 


Mq pc' 


-2 


1.85 


[2,1] 


Ensemble averaged clump density * 
CNM density * 




cm ^ 








?^CNM 


cm~^ 




32.9 


[1] 


WNM/ WIM density * 


l-WNM 


cm"^ 




0.349 


[1] 


CNM temperature * 


TcNM 


K 




85 


[1] 


WNM/ WIM temperature * 


TwNM 


K 




7860 


[1] 


FUV-flux ** 


X 


XD 








Metallicity ** 


z 


Zq 




1 




Clump mass spectral index ** 


a 






1.8 


[3] 


Mass-size relation index ** 


7 






2.3 


[3] 


Upper clump mass limit ** 




Mq 




10^ 


[3] 


Lower clump mass limit ** 


rni 


Mo 




10"^ 


[3] 



kpc), the molecular hydrogen gas mass surface density is 
assumed to have a constant value of 88.4 Mq pc~^ corre- 
sponding to an H2 mass of 10^ Mq within the central 600 
pc of the Galaxy (jPahmen et al.lll998l : iGiisten fc PhilippI 
|2004[) . This translates into a mean H2 gas mass density 
of 0.75 Mq pc^'^, corresponding to an H2 number density 
of about 15 cm^"^. The assumed mass distributions are vi- 
sualized in Fig. [2] a). The neutral atomic hydrogen mass 
distribution is used to estimate the contribution of diffuse 
atomic phases to the [C 11] 158//m fine structure line emis- 
sion (as shown in Section 14.2. ID. The contr ibution of the 
diffuse WIM is estimated bv lRevnoldj ()1990D to 37% of the 
total H I mass, which is slightly less than the WNM mass. 

We adopt a lower clump mass limit mi of 10^'^ Mq ac- 
cording to the observat ional lower clump mass limit from 
IHeithausen et all (Il998f). as well as an upper clu mp mass 



from IWolfire et al.l (|2003l ) (extrapolated inwards of 3 kpc 
Galactocentric radius) and the WIM d istributions are as- 
sumed to equal the WNM distributions (iMcKee fc Ostrikeil 
119771 : IRevnolds! [T991I : [Cro3l2005t iHiU et al.ll2007l) ! 



hmit m-u of 10^ Mq. iGorti fc HoUenbachl (|2002[ ) show an 



decreasing lifetime of clumps with decreasing mass and in- 
creasing density and argue against the existence of stable, 
small, and dense PDR clumps. This indicates a transient 
nature of these clumps. 

The observed mean gas density of molecular clouds 
shows variation w i th Galactocentric radius (e.g. 
iBrand fc WouterlootI I1995D . We adopt the Galactic 
radial distribution of the CNM density provided by 
IWolfire et cdl (|2003f) for the ensemble averaged clump 
density n^ns = Pons/mH , scaled to an absolute value of 
lO'^'^ cm"'^ at the solar circle (Fig. |2| b)), obtained by a 

fit to the intensity distributions as observed by COBE 
(see section 14. ip . As shown in Section I2.2[ the volume 
filling factor is determined by the underlying mass and 
density distributions. The result is shown in Fig. |2| c). 
The computed beam filling factors from the solar point of 
view along the Galactic plane are shown in Fig. |2] d) . The 
density and temperature distributions with Galactocentric 
radius of the CNM and WNM, used in the following to 
estimate the [C 11] emission from these phases, are adopted 



The FUV fiux in the Galaxy is mainly produced 
by OB stars. Observations show, that OB associa- 
tions are associated with gia nt molecular clouds (G MCs) 
(e.g. Stark fe Blitz 1978). McKee fc Williamsl (jl997t ) and 
IWilliams fc McKed (jl997D show, that the distributions of 
Galactic OB associations and GMCs can be described by 
truncated power l aws with the same sl ope, i.e. show a 
strong correlation. iMisiriotis et al.l ()2006D derive a nearly 
constant dust temperature from the COBE observations. 
Assuming a proportionality between dust continuum emis- 
sion tracing the molecular material and the incident FUV 
flux from the stellar population this leads to a rather con- 
stant average FUV-flux impinging on the molecular gas 
within the Galactic disk. According to lMcKee fc Williamsl 
(jl997l . and references therein) there are 6-7 OB associa- 
tions with more than 30 massive stars within 1 kpc dis- 
tance from the sun, dominati ng the local interstellar FUV 
radiation fleld. lBlaauwl (|l985l ) states a mean scale height of 
about 50 pc of the local OB associations. Therefore we as- 
sume one layer of surrounding OB associations with 3 next 
neighbors. This leads to an average radius of major influ- 
ence for each OB association of Rqb ~ 7~^/^ to 6~^/^ kpc 
or 378 to 408 pc. The actual distance measu rement to the 
Orion Nebula Cluster resulting in 414 ±8 pc l|Menten et al.l 
I2OO7D corresponds well to this ra nge. Togethe r with the 
local interstellar FUV-flux of 1 xu (|Drainelll978D and aver- 
age distances dMC-OB of about 20 to 50 pc of the GMCs to 
their associated OB associations (estima ted from the sizes 
and shapes given in iStark fc BlitafT97l Fig. 1 and Table 
1) we obtain FUV- fluxes impinging the molecular clouds of 
Xmc ~ (-RoB/rfMC-OB)^XD/3 with values of 19 to 139 xd 
or log(xMc/XD) ~ 1-3 to 2.1. In Section |3| we flnd a best 
fltting value of log(x/xD) ~ 1-8 for our model, falling well 
within this range. 



M. Cubick et al.: A clumpy-cloud PDR model of the global far-infrared line emission of the Milky Way 



7 




4.2 



4.0 



3.8 



o 3.6 



3.4 



-3.0 



-3.5 



> -4.0 



-4.5 



-5.0 




- 1.5 



180 




5 10 15 

Galactocentric Radius Rc [kpc] 




120 



60 -60 

Galactic longitude I [°] 



-120 



-180 



Fig. 2. a) Galactic mass surface density distributions of 
the molecular and atomic gas mass; b) ensemble averaged 
clump density distribution; c) resulting volume filling factor 
distribution with Galactocentric radius; d) derived beam 
filling factor distribution with Galactic longitude from the 
solar point of view. 




Fig. 3. Volume emissivity spectra of the clump ensemble 
for an ensemble averaged gas density of rigns = 10^ ® cm~'^ 
and an incident FUV-flux of x = 10^'^ xd- According to 
the logarithmic scaling of the clump-mass axis /(Md) = 
Mcid?7/dMci is plotted to see the integral contribution of 
each mass interval. The values are normalized to the func- 
tional value at Md = m\ for a convenient comparison of 
the different tracers. 



3. The clump-ensemble emissivity 

To discuss and understand the behavior of the clump en- 
semble emissivity it is useful to first look at the features 
of individual clumps and the mass-differential emissivity of 
the ensemble. As explained in Section [^751 the contribution 
to the volume emissivity spectrum of the ensemble drj/dMd 
in the various line transitions is dominated by the clump 
intensity Id, as the filling factors do not significantly de- 
pend on the clump mass. The logarithmic volume emissiv- 
ity spectrum d7^/dlog(A/ci) = Mcidrj/dMci for an ensem- 
ble averaged clump density of riens = W^'^ cm~'^ and an 
FUV-flux of X = 10^'^ xd is shown in Fig. [31 The symbols 
show the sampling points on the model grid. The upper 
plot shows the behavior for the different CO transitions. 
Note the trend towards larger contributions of the lower 
clump masses for higher transition numbers. The ensemble 
emission of the higher CO transitions (J > 5) almost exclu- 
sively stem from the low mass clumps. The bottom panel 
shows Mcid?7/dA/ci for the fine structure cooling lines. The 
[Cii] 158 iim and [Ci] 609 fim emissions show an increased 
contribution of high mass clumps. 
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Fig. 4. The CO emission of a single clump of one so- 
lar mass and a density of nd ~ lO'^'^ cm^'^ compared to 
the emission of an ensemble of clumps with the same total 
mass and ensemble averaged clump density, embedded in 
an isotropic FUV-flux of lO^ '^ xd- 

In Fig. [4] a direct comparison of the emission of an 
individual clump with the emission of a clump ensemble 
with 1 Mq is shown for the different CO transitions. For 
^3.8 



higher intensity in the mid- J CO transitions for the clump 
ensemble, which amounts to 2.5 orders of magnitude for 
the CO 8-7 line. This is due to the fact, that the mid-J 
CO emission shows a steep decrease with clump mass (see 
Fig. p. 

Fig.Oillustrates the dependence of the normalized cmis- 
sivity with variation of the lower and upper clump mass 
limits of the ensemble for different tracers. Again, the mid- 
J CO emission shows a strong dependence on the lower 
clump mass limit, corresponding to the strong intensity en- 
hancement with decreasing clump mass in these tracers. 
The CO 8-7 emissivity is nearly proportional to the lower 
clump mass limit. The [C ii] and [O i] emissivities show vari- 
ations of one order of magnitude over an upper clump mass 
limit range of three orders of magnitude. This is due to the 
increasing surface to volume ratio of the clump ensemble 
with decreasing upper clump mass limit (cf. Section 12. 2[) 
and the fact that the [C ii] and [O i] emissions arise mainly 
in shells in the outer parts of the clumps. 

4. Model results and comparison with observations 

The resulting emissivity distribution inherits the intensity 
dependence of the clumps on the Galactic mass and density 
distributions as shown in Fig. [51 

The beam averaged intensity is determined according 
to the derivation in Section [^751 For the comparison of our 
model result with the COBE FIRAS observations we cal- 
culate the intensities in a beam with an ex tent of 5° in 
Gala ctic longitude and 1° in Galactic latitude (jFixsen et al.l 
|1999[ ). All main PDR cooling hnes in the COBE FIRAS 
spectral range including [Cii] 158 //m, [Ci] 609 and 370 
/im, CO J = 1-0 to J = 8-7, and [O i] 146 /im are used for 
the comparison. 

To overview the emission features along the Galactic 
plane we labeled the main features visible in the [C ii] in- 
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Fig. 5. Dependencies of the clump ensemble emissivity rj 
on the upper and lower clump mass, mu respectively mi, 
for different tracers, Uens = lO"^'^ cm~^, and x = 10^'^ xd- 
The values are normalized to the upper mu = 10^ Mq 
respectively lower clump mass limit mi = 10^"^ -f^Qi 
used in the model. 



tensity cut in Fig. [Tj One can easily identify the Galactic 
large scale structure, dominated by the spiral arms and the 
Galactic molecular ring. 

4.1. Quantitative comparison 

To quantify the quality of the model result we derive the 
reduced x^ of the intensity distributions. The A/^ are the 
uncertainties of the observed intensities due to the cali- 
bration of the COBE FIRAS instrument and the data re- 
duction. The number of degrees of freedom is given by 
^dof = N\ Nspec — 2, with iVi the number of bins in Galactic 
longitude and Ngpec the number of species used for the fit. 
The subtraction of 2 is due to the 2 fitted parameters n^ns 
and X- Hence, the reduced x^ is given by 



1 



dof 



■ E 



ohs.i 



'mod, 2 



AL 



(24) 



The resulting x^-distribution of the intensity distributions 
is plotted as contours onto the density-FUV-plane in Fig. 
[8l The minimum x^ has a value of 16.5 at an ensemble aver- 
aged clump density at the solar circle of Ucns ~ 10^'^ cm~^ 
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Fig. 6. Calculated CO emissivity distribution vs. 
Calactocentric radius for ncns = lO'^'^ cm^'^ and x = 
10^'® XD at the solar circle. 
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Fig. 7. The [C ii] 158 /xm int ensity distribution with 
Calactic longitude as pubhshed bv lFixsen et al.l (|l999t ). The 
annotated source names are i d entifi ed using the CO 1-0 
map published bv lDame et al.l (|200ll ). 



and an FUV-flux of x = 10^'^ xd- The distribution ob- 
viously confines the density to FUV-flux ratio riens/x very 
well, as this ra tio determines the molecula r gas fraction 
within a PDR (jHoUenbach fc TielensI Il999t ). In contrast, 
the product nensX is hardly constrained. 

4.2. Discussion 

The results of the modeled line intensities are overlaid onto 
the observed intensities in Figs. [9l and [TOl We note that the 
relatively high value of the minimum is due to the pro- 
nounced Galactic large scale structure in the observed [C ii] 
emission (which is not contained in the smooth model dis- 
tribution) and the high signal to noise ratio of these data. 
In fact, it is surprising that one can find average values for 
these parameters, which allow us to reproduce the bulk of 
the FIR line emission of the Milky Way. The partly sig- 
nificant local deviations of the model results from the ob- 
servations are due to the strongly simplifying assumptions. 
Except CO J = 1-0, 6-5, and 8-7, which have large ob- 
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servational uncertainties, all lines are reproduced in total 
Galactic flux within factors below 2.3 (CO J = 2-1). The 
mismatch of the CO 2-1 indicates an additional cold molec- 
ular cloud component in the Galactic gas distribution. 

The resulting ^^CO 1-0 / CO 1-0 fine intensity ratio 
with a value of about 4 is in accordance with obse rvational 
results (iGierens et al.lll992l: iFaigarone et al.lll998[ ). 

4.2.1. Contributions to the Galactic [Cii] 158 /^tm emission 

Wc now turn to the comparison of the COBE observed [C ii] 
intensities with the model results (see Fig. Additional 
flux is measured mainly from the Galactic ring, the Cyg 
X region and other star-forming regions in the spiral arms. 
This is in contrast to the pictu re of [C ii]-emission stemmin g 
to a large part from the CNM (jHoUenbach fc Tiel ens"1999f). 
Based on dire ct UV absorption mea surem ents of the e xcited 
state of C+, iPottasch et all (|l979f ) and iGrv et al.l (|l992f) 
derive an average emissivity of the CNM, in particular dif- 
fuse Hi-clouds, of 5 X 10~^^ erg s~^ H-atom^^. Using this 
number the contribution from PDRs to the Galactic [C ii]- 
e mission should be very small. T his has been questioned 
bv lPetuchowski fc BennettI (|1993f ) who accounted for about 
half of the [C ii] -emission from PDRs, the other half from 
the WIM and only a v ery small fraction from the CNM. In 
contrast, iHeilej ()1994[) found that the CNM ranks second 
after the WIM wh ile PDR contributions should be small. 
ICarral et al.l (|1994f ) state that up to 30% o f the interstel- 
l ar [C ii] emission may stem from Hii- regions. iKramer et all 
(l2005t ) used ISO/LWS observations of the [Nil] 122 /xm line 
in the spiral arms of M83 and M51 to estimate that between 
15% and 30% of the observed [C ii] emission origin ates from 
Hii-regions. Computations by I Abel et al.l (|2005[ ) indicate 
that the main [C ii]-emission arises from PDRs and a frac- 
tion of 10% up to 60 % from H ii- regions. Calculations by 
iKaufman et"aLl (|2006t) find similar [C ii]-emission fractions, 
deviating by less than a factor of 2. 
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Fig. 10. The modeled intensities of the rotational transitions J = 1-0 (upper left) to J = 8-7 (lower right) of CO 
overlaid the observational distributions along the Galactic plane. Note that the combined CO J = 7-6 and [C i] 370 /im 
line intensities are shown in the lower left as they are not spectrally resolved by the COBE FIRAS instrument. 



We can test this result from the viewpoint of the ex- 
citation of C"*" which is a simple problem because the 
fine structure line arises basically from a two-level system. 
Using the excit ation rates for coll i sions with H, H2, He, and 
electrons from Flower fc Lau nayl (|1977[ ) , iLaunav fc Rouefj 
(|1977t) . andlwHionTBcli (.200l we can compute an upper 



limit to the cmissivity of C^ as a function of density under 
the temperature and ionization conditions of the different 
phases of the ISM via 



1 



nc+Au\Eu\(3 



An 



-EiiL 1 (I + IhlE) 



(25) 
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Fig. 9. The model results of the [Cii] 158 /im, the [Ci] 
609 /im, and the [O i] 146 /im intensity distributions with 
Galactic longitude overlaid on the observed distributions. 



with n the total number density of the collision partners 
(mainly hydrogen atoms or electrons), A^i the Einstein 
coefficient for spontaneous emission, (/i n the statistical 
weights for the lower and upper state, nQ+ the number 
density of E^i the energy spacing between the upper 
and lower state, n„ the critical density of the [C ii] ^P3/2" 
'^^1/2 transition with respect to neutral particles or ricr.e for 
electrons, ke the Boltzmann constant, and T the tempera- 
ture. The escape probability /3 describes the probabihty of 
a photon of the considered transition to escape the cloud 
from its position within, assuming constant excitation con- 
ditions allover the cloud. To obtain an upper limit of the 
emission of the atomic gas components, we assume opti- 



Table 5. The parameter values used for the escape prob- 
ability calculation of the [Cii] ^P3/2"^f 1/2 emissivity from 
the CNM and WNM. Powers of ten are given in parenthe- 
ses. 



Quantity 


Value 


/3 


1 


ffu 


4/3 


9i 


2/3 


A^i [s-i]^ 


2.4 (-6) 


Tier [cm"'^] 


4(3) 


Wcr.e [cm^^] 


4(1) 




Fig. 11. Optically thin emissivity of [Cii] under the ex- 
citation conditions of the different ISM components as a 
function of the gas density. 



cally thin emission, i.e. f3 ~ 1. The values for the different 
parameters arc summarized in Table [S] 

Fig. [11] shows the [Cii] emission from the different 
phases. All carbon is considered to be singly ionized, using 
an elemental abundance of 1.4 x 10~^ (see Table [3]). Upper 
limits are also guaranteed by assuming an ionization degree, 
i.e. an electron density relative to the total density of pro- 
tons, of 10-3 in the CNM, lO'^ in the WNM, and 1 in the 
WIM (see Table |4] for the local Galactic parameter values 
of the different phases). We find that the "standard value" 
of 5 X 10"^^ erg s~^ H-atom~^ from iHollenbach fc Tieleiii 
(I1999D is in agreement with the density range of the CNM 
given in IWolfire et al.l (j2003t ). At low densities an efficient 
excitation of C+ is only possible if a large fraction of the 
collision partners is charged, i.e. in the ionized medium. 

The results of iPetuchowski fc BennettI (|l993l ) andlAbel 
(|2006[ ) who found, that the [Cli] and [N 11] emission are 



strongly correlated, indicate that main contributions (60% 
maximum) are expected from low density Hii-regions. 
Using this result, we rerun the fit of the lines but increase 
the computed [C ii]-intensity by a global factor 2.5 before 
fitting the model results to the observations. This fit pro- 
vides a lower limit of the FUV-flux of x = 10^'^ xd, which 
is slightly below the estimated range from Section 12.41 with 
a lower limit of x = W^'^ Xd- The best fitting ensemble av- 
eraged clump density is derived to rions — lO"^'^ cm"^ in this 
case. Additionally we check the fit results completely ignor- 
ing the [C 11] intensity distribution. In this case we obtain 
parameter values of no„s = lO^'^ cm^'^ and x = lO"'^" xd- 
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Fig. 12. PDR [Cii] intensity corresponding to the best fit 
parameter values ignoring the [C ii] emission and under the 
assumption that the contribution of the PDR-correlated 
Hii regions to the total [Cii] emission is 60%. 
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Fig. 13. The modeled FIR intensity distribution along 
the Galactic plane in the COBE FIRAS wavelength band 
overlaid on the o bserved distribution as published by 
iFixsen et all (|l999t) . 



5. Summary and Conclusions 



The x^-distribution with a minimum of 5.40 for this case 
is also shown in Fig. [8] The resulting low FUV-flux is due 
to the low signal to noise ratio in the [O i] and upper CO 
transitions, so that this result is biased by the low-J CO 
and [C i] emission and lies below our estimate for the av- 
erage FUV-flux seen by the Galactic molecular gas. The 
confinement of the distribution ignoring the [C ll] inten- 
sity distribution is weak in cither direction. In particular 
at rions = 10^'^ cm~^ and x ~ 10^'^ xd we find a value 
of 5.75 - only about 7% more than the minimum value, 
indicating that the [C ii] emission fixes the resulting FUV- 
flux value of the model. The computed PDR [C ii] intensity 
distributions for these cases are shown in Fig. [T^ 

There are smaller, comparable contributions from the 
diffuse WIM, which has an average density well below 1 
cm-3 (iMcKee fc Ostrikeil ll977l: iRevnoldd [Tool ICro3l2005t 



iHill et al.ll2007 ). and from the CNM , and a vanishing con- 
tribution of the WNM. However, the much larger scale 
height of the diffuse WIM leads to a reduced contribution 
in comparison to the CNM, in agreement with our result for 
the intensity distribution along the Galactic plane shown in 
Fig. [3 

We conclude that our model can reproduce the Galactic 
[C ii] emission with the main contribution from clumpy 
PDRs, but we cannot exclude a major contribution from 
low-density Hii regions. 



The bulk of the FIR and submillimeter line emission in 
[Cii] 158 nm, [Ci] 609 fim and 370Atm, [Oi] 146 fim, and 
CO 1-0 to 8-7 of the Milky Way as observed by COBE 
FIRAS can be reproduced by a clumpy PDR-model, which 
takes the fractal structure of the dense ISM as a clump 
ensemble with a power law clump mass and size distribu- 
tion into account. This is a remarkable result, as the model 
contains essentially no free parameters; all parameters be- 
ing constrained by independent knowledge about the struc- 
ture, kinematics of, and star formation rate in the disk of 
the Milky Way. 

The result suggests that the bulk of the Galactic FIR 
line emission stems from fractally structured PDRs, illumi- 
nated by the FUV radiation field emitted by the Galactic 
population of massive stars in adjacent star forming regions. 

The UV-penetrated, clump-cloud PDR scenario applied 
here to the Milky Way line emission as observed by COBE 
is a versatile approach, which can be applied in the future 
to model the emission from individual Galactic star forming 
regions as well as the line maps of spatially resolved nearby 
galaxies. 

Future studies and observations should address the va- 
lidity of particular assumptions, like the limits of the clump 
mass distribution. Observations with high spatial and/or 
spectral resolution, that will be possible with upcoming 
missions as SOFIA, Herschel, and ALMA, will give further 
insights in the nature of the ISM. 



4.2.2. The FIR continuum emission 

Only 51% of the Galactic FIR emission in the COBE 
FIRAS spectral range is reproduced by a separate dust 
radiative transfer calculation for the KOSMA-r PDR- 
model clumps (R. Szczerba, priv. comm.) as shown in Fig. 
fT3l This is expected, as the dust is also heated by less ener- 
getic radiation than the FUV f rom t he older stellar popula- 
tion. iMochizuki fc Nakagawal (j200GD argue that about half 
of the total stellar flux lies in the FUV spectral range in 
agreement with our result. 
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